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ABSTRACT 

Administration of mesenchymal stem cells (MSCs) 
has the potential to ameliorate degenerative dis- 
orders and to repair damaged tissues. The homing 
of transplanted MSCs to injured sites is a critical 
property of engraftment. Our aim was to identify 
microRNAs involved in controlling MSC proliferation 
and migration. MSCs can be isolated from bone 
marrow and umbilical cord Wharton's jelly (BM- 
MSCs and WJ-MSCs, respectively), and WJ-MSCs 
show poorer motility yet have a better amplification 
rate compared with BM-MSCs. Small RNA seq- 
uencing revealed that miR-146a-5p is significantly 
overexpressed and has high abundance in WJ- 
MSCs. Knockdown of miR-146a-5p in WJ-MSCs 
inhibited their proliferation yet enhanced their mi- 
gration, whereas overexpression of miR-146a-5p 
in BM-MSCs did not influence their osteogenic 
and adipogenic potentials. Chemokine (C-X-C 
motif) ligand 12 (CXCL12), together with SIKE1, 
which is an l-kappa-B kinase epsilon (IKK/) sup- 
pressor, is a direct target of miR-146a-5p in MSCs. 
Knockdown of miR-146a-5p resulted in the down- 
regulation of nuclear factor kappa-B (NF-kB) 
activity, which is highly activated in WJ-MSCs and 
is known to activate miR-146a-5p promoter. miR- 
146a-5p is also downstream of CXCL12, and a 



negative feedback loop is therefore formed in 
MSCs. These findings suggest that miR-146a-5p is 
critical to the uncoupling of motility and proliferation 
of MSCs. Our miRNome data also provide a 
roadmap for further understanding MSC biology. 



INTRODUCTION 

Human mesenchymal stem cells (MSCs) have been 
identified as multipotent mesoderm-derived stromal cells 
that have the ability to self- renew and differentiate (1); 
they have been applied as clinical treatments for bone 
and other tissue defects (2-A). On activation by tissue 
damage, MSCs contribute to tissue-repair processes 
through a multitude of activities, including cell prolifer- 
ation, migration and differentiation. The mobilization of 
bone marrow (BM)-derived MSCs from BM to the periph- 
eral blood, and their eventual entry into the injured brain, 
plays a crucial step in brain plasticity and stroke therapy 
(5). MSC activities also affect the therapeutic efficacy of 
engraftment, specifically if only low number of transplanted 
MSCs migrate to the injured site after infusion, which will 
limit the therapeutic applications of MSCs (6). The expan- 
sion/proliferation rate of MSCs ex vivo also influences cell 
motility, as MSCs lose their mobility during cultivation (7). 

microRNAs (miRNAs) are short non-coding RNAs 
(~22nt) that can repress translations through imperfectly 
binding to target messenger RNA. After being transcribed 
and processed by Drosha and Dicer, miRNAs are then 
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loaded into an RNA-induced silencing complex that brings 
about the regulation of translation (8). To date, relatively 
few studies have examined miRNA functionality in MSCs: 
miR-335 has been shown to inhibit cell proliferation, mi- 
gration and differentiation (9). In addition, miR-138 modu- 
lates osteogenesis by MSCs (10). miR-204 has also been 
found to inhibit osteogenesis but to promote adipogenesis 
by MSCs (11). We recently found that miR-34a is able to 
modulate the cellular motility genes of neural precursor 
cells derived from Wharton's jelly MSCs (WJ-MSCs) (12). 

Till date, hundreds to thousands of miRNAs have been 
identified in animals and plants, and many more miRNAs 
are continuously being identified by newly available 
technologies, including small RNA sequencing (smRNA- 
Seq). High-throughput sequencing is able to not only reveal 
the expression profiles of known miRNAs but also identify 
other non-coding small RNAs and discover new miRNAs 
that have not been recorded previously in any databases, in 
particular the miRBase repository. smRNA-Seq has been 
used to carry out research on various types of stem cells, 
including embryonic stem cells (13-15), hematopoietic stem 
cells (16) and neural precursor cells (13). Novel miRNAs 
have also been identified using smRNA-Seq during neural 
differentiation of embryonic stem cells (15) and during 
endothelial differentiation (17). Nevertheless, no 
smRNA-Seq work has been reported on somatic MSCs. 

Because the implanted number and homing of trans- 
planted MSCs to injured sites is one of the critical properties 
in relation to engraftment, in the present study our aim was 
to identify miRNAs that are involved in controlling the 
proliferation and migration phenotypes of MSCs. 
We hypothesized that miRNAs involved in stem cell 
motility and proliferation must be present in undifferenti- 
ated MSCs, given the variations observed on their 
mobility. MSCs from different sources have different func- 
tionality. MSCs can be obtained from BM as well as other 
fetal or postnatal tissues, including adipose tissue, umbilical 
cord blood and the Wharton's jelly of the umbilical cord 
(18). WJ-MSCs have been considered as a great alternative 
source for the harvesting of MSCs ( 1 9) and have multilineage 
differentiation ability that allows them to become osteocytes, 
chondrocytes, adipocytes, cardiomyocytes and neurons (20). 
The differentiation, immunomodulation and proliferation 
abilities, as well as transcriptome profiles, of WJ-MSCs 
have been compared with other MSCs (21-23). BM-MSCs 
have better osteogenic and adipogenic abilities, whereas WJ- 
MSCs have a higher proliferation potential (21). The migra- 
tion capacity of BM-MSCs was found higher than that of 
WJ-MSCs (24). In the present study, we applied smRNA- 
Seq to WJ-MSCs and BM-MSCs for identifying miRNAs 
that might regulate the activation of repair-related pheno- 
types of MSCs, as well as various specific target genes 
associated with these phenotypes. 

MATERIALS AND METHODS 

Isolation and cultivation of MSCs from BM and the 
Wharton's jelly of umbilical cord 

All human BM and WJ cells were obtained from unrelated 
donors with informed consent, and this study follows the 



guidelines of 1975 Helsinki Declaration. BM-MSCs and 
WJ-MSCs were isolated and cultured as previously 
described (21). In brief, all MSCs were cultured in 
MesenCult (Stem Cells Technology) or MesenPRO RS 
(Invitrogen) medium at 37° C in a 5% C0 2 humidified in- 
cubator. Mononuclear cells from BM were separated by 
centrifugation using a Histopaque-1077 gradient (1077-1; 
Sigma) and plated in fresh culture medium. WJ cells were 
isolated from the umbilical cord using collagenase and 
trypsin within 24 h of collection. In both cases, non- 
adherent cells were removed after 24-72 h plating. On 
reaching confluency, the cells were treated with 0.05% 
trypsin-ethylenediaminetetraacetic acid and subcultured 
in fresh culture medium. 

smRNA-Seq and data analysis 

smRNA-Seq was performed on the Illumina Solexa GA 
IIx sequencing system (Illumina, San Diego, CA, USA) 
according to manufacturer's instruction. The results 
were analyzed by an in-house bioinformatics pipelines 
(25). In brief, according to the length and sequence, iden- 
tical adapter-trimmed Fastq sequences were integrated to 
give a series of non-redundant sequences and aligned to 
genome [University of California, Santa Cruz (UCSC) 
hgl9] using the Bowtie aligner (26). Sequences that were 
able to be aligned with the transcribed genome region of 
known miRNAs (according to miRBase R19, http://www. 
mirbase.org/) were quantified as absolute expression 
values for each known miRNA. Expression values of 
specific miRNAs were calculated as the RPM value 
(reads per million mapped reads), namely RPM = C/ 
MN x 10 6 , where C is 'read numbers aligned to a given 
miRNA chromosomal region', M is 'multiple mapping 
numbers across all miRNA regions of given read (i.e. 
the number of different miRNA chromosome locations 
mapped by the same read)' and N is 'total read numbers 
that map to human genome sequence of a single sample'. 
To minimize the effect of cross-mapping of sequences with 
uncertained genomic locations, the expression is divided 
by its number of cross-mapping events. The q-value (false 
discovery rate-adjusted P-value) was done by 
edgeR package (27) of the Bioconductor (http://www. 
bioconductor.org) suite of software for the R statistical 
programming language (http://www.r-project.org). 

All smRNA-Seq data were deposited in NCBI Gene 
Expression Omnibus database (accession number 
GSE46989). 

Microarray analysis 

The total RNA from BM-MSCs transduced with pCDH- 
146a-5p or pCDH vector, as well as those of WJ-MSCs 
transfected with short interfering RNAs targeting miR- 
146a-5p or a green fluorescent protein (GFP) control, 
was isolated. mRNA profiling was carried out using 
Affymetrix Human U133 Plus 2.0 chips, and data 
analysis was performed as described previously (28). 
Robust Multi-array Average (RMA) log expression units 
were calculated from Affymetrix GeneChip array data 
using the 'affy' Bioconductor package in the R 
programming language. Gene expression levels were 
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normalized using default RMA settings. The array results 
are deposited in the Gene Expression Omnibus database 
with accession number GSE44186. 

Quantification of mRNA and miRNA expression 

RNA extraction and real-time quantitative polymerase 
chain reaction (RT-qPCR) were performed as described 
previously (12). Total RNA was extracted using Tri 
Reagent (Sigma-Aldrich Co., St. Louis, USA) according 
to the manufacturer's instructions. For mRNA qPCR, 
complementary DNA was synthesized from 100-1000 ng 
of total RNA using RevertAid™ Reverse transcriptase 
kit (K1622, Fermentas, MA, USA). Amplification of 
genes was performed using specific primers 
(Supplementary Table S2), and Maxima™ SYBR Green 
qPCR Master Mix (K0222, Fermentas), on the 
StepOne™ sequence detector (Applied Biosystems, 
USA). The relative expression levels of target genes were 
normalized against the expression level of glyceraldehyde 
3-phosphate dehydrogenase (GAPDH). For miRNA 
qPCR, the expression levels of specific miRNAs were 
detected using stem-loop RT-qPCR (29). The primers 
designed and used for each miRNAs are listed in 
Supplementary Table S2; snRNA U6 was used as an 
internal control for normalization of miRNA expression. 

Dual luciferase reporter assays 

The full length of the 3' untranslated region (3' UTR) of 
SIKE1, and a 2-kb fragment from the 3' UTR of CXCL12 
that contains the predicted target site, as well as their 
various mutant forms that had been mutated within 
their seed regions, were cloned into the pGL3 basic 
reporter plasmid (Promega, Charbonnieres, France). 
Next, 293T cells were co-transfected with the firefly 
luciferase reporter construct (500 ng) and either the miR- 
146a expression plasmid or the vector control (500 ng), 
together with an internal Renilla pRL-TK luciferase 
control plasmid (50 ng). The transformed cells were har- 
vested at 24 h after transfection and their luciferase 
activity measured as the ratio of firefly luciferase activity 
to Renilla luciferase activity. 

MSC migration, proliferation and differentiation assays 
Cell migration 

A total of 5 x 10 4 MSCs in 100 ul of culture medium were 
seeded into the upper insert with an 8-um pore-size 
membrane (Corning, NY, USA), whereas 600 ul fresh 
culture medium containing 10% serum was added to the 
lower chamber. After 12 h of incubation, each membrane 
was fixed with 4% paraformaldehyde, and the MSCs 
present on the membrane were then stained with 
Hoechst 33342 solution (Sigma). The various numbers of 
migrated MSCs were determined by counting the number 
of nuclei beneath the filter membrane. 

Cell proliferation 

BM-MSCs and WJ-MSCs had their incubation started at 
the same passage, and the number of cells was counted 
every five days after subculture. For the 3-(4,5- 
dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide 



(MTT) proliferation assays, WJ-MSCs were transduced 
with the appropriate lentivirus expression vector. The 
day of infection was set as day 0. After various times of 
incubation, the cells were treated with 1 % thiazolyl blue 
tetrazolium for 1-2 h at 37°C, followed by 0.1% sodium 
dodecyl sulfate in 2-propanol and thorough mixing. The 
results were obtained by measuring absorbance at a 
wavelength of 570 nm using a multiwell scanning spec- 
trophotometer (Thermo, Multiskan Spectrum). 

Cell differentiation 

MSC differentiation was performed as described previ- 
ously, with some modifications (21). BM-MSCs were 
transduced with lentivirus expressing miR-146a-5p or the 
vector control. One day post-infection, the medium was 
replaced with osteogenic or adipogenic induction medium. 

RNA interference and plasmids of miRNA overexpression 
and knockdown 

miR-146a-5p antagomir (IH-300630-05-0005; Dharma- 
con, MA, USA) or control oligomers (siRNA against 
GFP, Sigma, 20 nM) were transfected using a Micro- 
porator MP- 100 (Invitrogen, USA). miR-146a-5p expres- 
sion plasmid was purchased from SBI (System 
Biosciences, CA, USA). The lentiviral expression vector 
expressing hairpin sequence against miR-146a-5p was 
also purchased from SBI (System Biosciences, miRZip, 
MZIP000AA-1). 

Flow cytometry and immunoblotting 

Flow cytometry and immunoblotting were performed as 
described previously (21). For flow cytometry, MSCs in 
phosphate buffered saline were incubated with mono- 
clonal antibodies against CD29 (303004, Biolegend), 
CD34 (MCA547PE, Serotec), CD44 (312306, Biolegend), 
CD45 (304006, Biolegend), CD73 (550257, BD), CD90 
(MCA90F, Serotec) or HLA-DR (307604, Biolegend). 
Unstained MSCs were used as a negative control. For 
immunoblotting, antibodies against phospho-S536 p65 
(#3031, Cell Signaling, USA), p65 (06-418, Millipore, 
USA) or GAPDH (Millipore, USA) were used, and 
GAPDH acted as the internal control. 

Statistical analysis 

All results are presented as mean ± standard deviation 
with statistical significance tested using 2-tailed Student's 
z-tests. 

RESULTS 

Global analysis of miRNA expression in human MSCs by 
sequencing small RNAs 

We isolated human MSCs from BM and WJ as described 
(21). Both types of MSC had a similar morphology and 
expressed common surface markers, namely CD73, 
CD29, CD90 and CD44 (Supplementary Figure SI). 
Hematopoietic markers such as CD34 and CD45 were 
lacking in both populations (Supplementary Figure SIB). 
We have previously shown that BM-MSCs show better 
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Figure 1. Differential proliferative and migratory abilities between 
BM-MSCs and WJ-MSCs. (A) Proliferative potentials of two MSCs. 
The total cell numbers of cultured WJ-MSCs and BM-MSCs were 
counted every 5 days till day 20 after plating. Results are presented 
as mean ± SD in two duplicates. *P<0.05 by (-test. (B) Migratory 
ability of WJ-MSCs and BM-MSCs. Migrated cells were counted 
after 12 h incubation. Fold change in migration ability is compared 
against the WJ-MSCs. Representative images of migrated MSCs after 
staining with Hoechst 33342 are shown. Results are presented as 
mean ± SD from two independent experiments. 



osteogenic/adipogenic differentiation abilities than WJ- 
MSCs (21). We further evaluated functional variations in 
relation to the BM-MSCs and WJ-MSCs by subjecting cells 
of the same passage to cell proliferation and cell migration 
assays: WJ-MSCs were found to proliferate faster, whereas 
BM-MSCs had better cell motility (Figure 1). 

We hypothesized that stem cell miRNAs are involved in 
regulating these phenotypes. To explore this, we per- 
formed smRNA-Seq to decipher the expression levels of 
non-coding small RNA, especially miRNAs, expressed by 
BM-MSCs and WJ-MSCs; these MSCs had been obtained 
from four and three different donors, respectively, for 
allowing us to ignore inter-individual variations. To 
filter out sternness miRNAs that might be involved in 
MSC motility and/or growth, we also compared the 
miRNome patterns of MSCs with those of differentiated 
osteocytes and adipocytes derived from BM-MSCs. The 
successful differentiation of BM-MSCs was confirmed by 
the analysis of various markers by qPCR (Supplementary 
Figure S2). Using the Illumina Solexa platform, each 
sample smRNA-Seq reaction generated 17-36 million of 
high-quality reads that corresponded to 87-586 thousand 
of non-redundant reads (Supplementary Table SI). 

Patterns of known miRNAs in undifferentiated and 
differentiated cells 

A total of 12 stem cell miRNAs and 23 progeny miRNAs 
were found to be differentially expressed when undifferen- 
tiated and differentiated cells were compared. A heatmap 
was then created to show details of their differential 




B 



2.0 



c 
o 



in 
cd 



O 



o 



>. 
o 
o 

CD 



o 

Q_ 
< 



* i.o 

> 

1 0.5 

CD 

a. 

0.0 



miR-503 



WJ BM 



c 
o 
'35 

"> . 
2 4 

Q. 
X 
II) 

CD 

> 2 
+j 

CD 

a. 



o 

miR-345-5p 



JZL 



n 



2.5-. miR-423-5p 



WJ BM 



c 

I 2 

(A 
CD 

Q. 1.5- 

X 
CD 

a 1.0- 

1 o 

0.0 



a 



WJ BM 



Figure 2. miRNA patterns in undifferentiated MSCs and differentiated 
cells derived from BM-MSCs. (A) A heatmap shows the differential 
expression of miRNAs among stem cells and differentiated progeny 
cells. miRNAs differentially expressed between undifferentiated and 
differentiated cells (<j<0.01) were shown. (B and C) Validation of dif- 
ferential expressed miRNAs by RT-qPCR. Highly expressed miRNAs 
in undifferentiated MSCs (B) or differentiated cells (C) are shown. O: 
osteocytes; A: adipocytes. 



expression across stem cells and differentiated cells 
(Figure 2A), and these findings are expanded in 
Supplementary Figure S3. Examples of the miRNAs that 
were found to be upregulated in stem cells (miR-503) or in 
osteocytes/adipocytes (miR-345-5p, miR-423-5p) were 
validated by RT-qPCR (Figure 2B and C). 

A total of 68 cell type-specific known miRNAs were 
identified as having fold changes of >2 and RPM >20: 6 
miRNAs were found to be unique to BM-MSCs, 
40 miRNAs to WJ-MSCs, 15 to osteocytes and 7 to 
adipocytes. A heatmap and details of their differential 
expression are presented in Figure 3A. Among all cell 
type-specific miRNAs in BM-MSC and WJ-MSC, miR- 
146a-5p was identified as significantly overexpressed at 
high abundance in WJ-MSCs (Supplementary Figure 
S4). RT-qPCR verified the unique expression patterns of 
miR-146a-5p and miR-29b-3p in WJ-MSCs, and miR- 
140-3p and miR-193-5p in BM-MSCs (Figure 3B). The 
expression level of miR-146a-5p was not significantly 
changed during MSC osteogenic or adipogenic differenti- 
ation (Supplementary Figure S5). 

Functional module analysis reveals a crucial role for miR- 
146a-5p in MSC migration and proliferation 

The above profiling data imply that miR-146a-5p is a stem 
cell-specific miRNA that is likely to be involved in MSC 
motility and/or proliferation, but not in differentiation. 
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Figure 3. Cell type-specific miRNAs. (A) A heatmap shows the specific expression of miRNAs in each cell type. (B) Validation of miRNA levels in 
different MSCs by RT-qPCR. Each spot represents an independent donor. Data are shown as mean of at least four independent donors. *P<0.05 
by (-test. 



We used the Ingenuity Pathway Analysis (IPA) software 
to evaluate the possible functions of miR-146a-5p in 
MSCs. A total of 334 BM-MSC-enriched genes (from 
our previous gene expression array study (21)) were sub- 
jected to IPA to pinpoint possible miR-146a-5p function- 
ality in MSCs. Genes regulated by miR-146a-5p 
(discovered by either wet lab experiments or bioinfor- 
matics target prediction) were connected to miR-146a-5p 
by IPA (Supplementary Figure S6). Genetic network 
analysis showed that miR-146a-5p targets genes were 
involved in both cell proliferation and migration, 
with P-values of 5.69E-3 and 4.40E-3, respectively 
(Supplementary Figure S6). 

miR-146a-5p (previously known as miR-146a) has been 
associated with nuclear factor kappa-B (NF-kB) activity 
(Supplementary Figure S6); specifically, expression of 
miR-146a-5p is regulated by NF-kB through NF-kB 
binding sites in miR-146a-5p promoter (30), whereas NF- 
kB is targeted by miR-146a-5p in breast cancer (31); these 
interactions form a negative feedback loop. miR-146a-5p is 
also crucial for stem cell properties whereby miR- 146a-5p is 



able to directly induce hematopoietic stem cells to differen- 
tiate into macrophages, as well as being able to modulate 
neural stem cell differentiation (32,33). miR-146a-5p has 
also been shown to induce the proliferation of vascular 
smooth muscle cells (34). Furthermore, miR-146a-5p is 
able to suppress cell invasion in pancreatic cancer (35). 
However, among MSCs, the various roles of miR- 146a- 
5p in NF-kB regulation, cell proliferation and migration 
have not been investigated. 

miR-146a-5p regulates MSC migration and proliferation, 
but not differentiation 

Because miR-146a-5p was found to be more abundant in 
WJ-MSCs than in BM-MSCs, we analyzed the effects of 
miR-146a-5p on MSC migration and proliferation. WJ- 
MSCs transfected with miR-146a-5p antagomirs showed 
a 50% reduction in expression level of miR-146a-5p, 
whereas in parallel their cell motility was increased 
compared with cells transfected with the control 
oligomer (si-Control in Figure 4A). BM-MSCs transduced 
with lentiviruses expressing miR-146a-5p at a multiplicity 
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Figure 4. miR-146a-5p regulates MSC migration and proliferation. (A 
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miR-146a-5p (shl46a), and cell proliferation rate was measured at the 
indicated time points using MTT assays. Results are presented as 
mean ± SD in two duplicates. *P<0.05 by Mest. 



of infection of five showed an increased level of miR-146a- 
5p, and this was accompanied by decreased cell motility 
(Figure 4B). These results suggest that miR-146a-5p ex- 
pression controls MSC motility. 

To determine the role of miR-146a-5p in MSC prolifer- 
ation, MTT assays were conducted on WJ-MSCs at 
various time points after miR-146a-5p was knocked 
down. Cell proliferation was reduced in WJ-MSCs 
transduced with a lentivirus expressing a hairpin 
sequence targeting miR-146a-5p (shl46a) (Figure 4C). 
Previous studies have shown that NF-kB plays a crucial 
role in regulating MSC survival, proliferation and migra- 
tion: Toll-like receptor ligand is able to activate NF-kB 
activity, reduce MSC migration ability and increase MSC 
proliferation, which implies that the NF-kB pathway con- 
tributes in MSC functionality (36). We have previously 



shown that the NF-kB protein level is higher in WJ- 
MSCs than in BM-MSCs (21). We found that knockdown 
of miR-146a-5p in WJ-MSCs led to a significant reduction 
in p65 phosphorylation (Supplementary Figure S7A) 
when examining NF-kB activity in WJ/shl46a or WJ/ 
shControl cells using immunoblotting with antibodies 
against the phosphorylated/activated p65 (RelA) subunit 
of the NF-kB heterodimer. These results suggest that an 
increase in miR-146a-5p level in MSCs is able to stimulate 
NF-kB activity. 

The various roles of miR-146a-5p in MSC differenti- 
ation were also monitored. BM-MSCs transduced with 
miR-146a-5p lentivirus or vector control were subjected 
into osteogenic and adipogenic differentiation, and 
various lineage-specific differentiation markers were 
measured every week up to week 2 or 3 using RT-qPCR 
or reverse transcription polymerase chain reaction 
(RT-PCR). The expression of differentiation markers 
was found to be unaffected by miR-146a-5p overexp- 
ression (Supplementary Figure S7B and C), which 
implies that miR-146a-5p does not participate in MSC 
osteogenesis or adipogenesis. 

miR-146a-5p targets CXCL12 and SIKE1 in MSCs 

We used a combination of bioinformatics target predic- 
tion and gene-expression microarray data to pinpoint 
potential miR-146a-5p downstream targets to gain 
mechanistic insights into miR-146a-5p functionality. The 
transcription profile of BM-MSCs with miR-146a-5p 
overexpression, as well as that of WJ-MSCs with miR- 
146a-5p knocked down, was explored by gene-expression 
microarray. We focused on genes that were downregulated 
in BM-MSCs when miR-146a-5p was overexpressed, and 
were upregulated in WJ-MSCs when miR-146a-5p was 
knocked down (fold change >1.5). A total of 73 genes 
were identified at this step. Next, because miR-146a-5p 
has been shown to be downregulated in BM-MSCs, we 
considered only genes that were more abundant in BM- 
MSCs. A total of 15 of the original 73 genes fitted this 
criterion. Among these 15 genes, 7 possessed putative 
miR-146a-5p binding sites in their 3' UTRs according to 
the miRanda and PITA prediction programs. Among 
these genes, we focused on CXCL12 and SIKE1; this 
was because CXCL12 is able to induce MSC migration 
(37), and SIKE1 is a suppressor of IKKs, a NF-kB modu- 
lator through a non-canonical pathway (38). We found 
that the endogenous levels of CXCL12 and SIKE1 were 
higher in BM-MSCs than in WJ-MSCs (Figure 5A), and 
that the mRNA expression levels of CXCL12 and SIKE1 
were increased in WJ-MSCs that had been transfected 
with miR-146a-5p antagomirs using cells from independ- 
ent individuals (Figure 5B). 

The direct repression of CXCL12 and SIKE1 by miR- 
146a-5p was proved by reporter assays. Specifically, miR- 
146a-5p was found to repress expression of luciferase in 
constructs containing the CXCL12 or SIKE1 3' UTR 
fused downstream of the luciferase gene, respectively 
(Figure 5C and D, the wild type (WT) group). Thus, the 
luciferase assays showed that miR-146a-5p directly inter- 
acts with the CXCL12 and SIKE1 3' UTRs, which results 
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Figure 5. CXCL12 and SIKE1 are direct targets of miR-146a-5p. (A) The relative expression levels of CXCL12 and SIKE1 in both MSC types were 
detected by RT-qPCR (n = 4). (B) CXCL12 and SIKE1 mRNA levels in WJ-MSC transfected with miR-146a-5p antagomirs were analyzed by RT- 
qPCR. Results are shown as mean ± SD from two duplicates. *P<0.05 by Mest. (C and D) Putative miR-146a-5p target sites in the CXCL12-3' 
UTR (C) and SIKE1-3' UTR (D). UTR: untranslated region; WT: wild type; MUT: mutant. (E and F) Luciferase reporter assays using the CXCL12- 
3' UTR (E) or the SIKE1-3' UTR (F) in 293T cells transfected with control (Vec) or miR-146a-5p expression vector. Results are presented as 
mean ± SD of three independent experiments. 



ill a modification of gene expression (Figure 5E and F). 
This repression was reversed by independently mutating 
the putative miR-146a-5p-binding sites in the CXCL12 
and the SIKE1 3' UTR (Figure 5C and D, the Mut 
group). In the SIKE1 3' UTR, only mutation of sitel 
restored luciferase expression, showing that sitel is the 
potential binding site for miR-146a-5p (Figure 5F). 

miR-146a-5p modulates MSC motility and proliferation 
via CXCL12 and SIKE1 

To further dissect the functional impact of CXCL12 and 
SIKE1 on miR-146a-5p-mediated MSC motility and pro- 
liferation, we compared the migration and amplification 
abilities of MSCs that had either CXCL12 or SIKE1 
knockdown. Knockdown of either CXCL12 or SIKE1 
was able to significantly decrease the migration ability of 
WJ-MSCs (Figure 6A and B). The proliferation rate of 
WJ-MSCs was reduced when the miR-146a-5p level was 
knocked down. Furthermore, knockdown of SIKE1 in the 
miR-146a-5p knockdown WJ-MSCs, but not knockdown 
of CXCL12, was able to restore WJ-MSC proliferation 
(Supplementary Figure S8A). 

We next examined whether the proliferation rate of 
WJ-MSCs was associated with the phosphorylation level 
of the NF-kB p65 subunit. p65 phosphorylation was 
found to be decreased in miR-146a-5p-knockdown 
WJ-MSCs, but this reduction could be rescued in MSCs 
when there is miR-146a-5p and SIKE1 double knockdown 
(Supplementary Figure S8B). Thus, miR-146a-5p 



enhanced NF-kB activity in WJ-MSCs and, furthermore, 
NF-kB signaling is able to induce miR-146a-5p expression 
through binding to the NF-kB binding sites that are 
present in the miR-146a-5p promoter (30). Knockdown 
of CXCL12 resulted in an increased expression of miR- 
146a-5p in WJ-MSCs (Figure 6C), implying that miR- 
146a-5p positively regulated its expression in MSCs via 
CXCL12 and NF-kB signaling (Figure 7). 



DISCUSSION 

Cell-based therapies using multipotent MSCs for organ 
regeneration have been and are being pursued for a 
variety of diseases. Clues on molecules controlling MSC 
biology can be obtained by comparing MSCs with differ- 
ent biological activities. For example, a comparison of 
MSCs isolated from Murphy Roths Large MRL/MpJ 
mice was able to demonstrate enhanced regenerative 
capacity compared with MSCs from C57BL/6 (WT) 
mice; these results suggest that sFRP2 is a key molecule 
that affects biogenesis and is important to obtaining a 
superior MSC regenerative phenotype (39). It is clear 
that MSCs from different human anatomic locations 
possess diverse biological activities; these include differ- 
ences in differentiation, proliferation and migration po- 
tentials. miRNA functions are suggested to regulate 
various biological processes in stem cells (40,41). To thor- 
oughly characterize the miRNA changes in MSC, in this 
study we used smRNA-Seq to profile miRNA expression 
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Figure 6. Knocking down CXCL12 and SIKE1 inhibits WJ-MSC 
migratory ability, and CXCL12 silencing increases miR-146a-5p 
levels. (A and B) CXCL12 and SIKE1 control MSC motility. A total 
of 5 x 10 4 WJ-MSCs were transduced with CXCL12 (A) or SIKE1 (B) 
shRNA and then assessed by Transwell migration assay. Migrated cells 
were counted after 12 h. Results are shown as mean ± SD from two 
duplicates. */ > <0.05 by Mest. (C) miR-146a-5p expression in WJ- 
MSCs after knockdown of CXCL12. WJ-MSCs transduced with 
CXCL12 shRNA for 48 h were harvested. All results are presented as 
mean ± SD from two duplicates experiments. *P < 0.05 by /-test. 
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patterns in WJ-MSCs, BM-MSCs and their differentiated 
progenies. Follow-up wet lab studies showed that MSCs 
use a complex miRNA-mRNA network to modulate their 
functionalities. In this network, miR-146a-5p and its 



downstream genes play a crucial role in inhibiting cell mi- 
gration while promoting cell proliferation. 

Previous studies have demonstrated that miRNAs act as 
important regulators for MSCs. For example, the expres- 
sion of miR-193 was increased in low-level laser irradi- 
ation-treated MSCs, which enhanced MSC proliferation 
by targeting inhibitor of growth family, member 5 (42). 
Gain-loss function experiments showed that miR-21 nega- 
tively regulates the proliferation of adipose tissue-derived 
MSCs through regulating signal transducer and activator 
of transcription 3 (STAT3) expression (43). As to cell dif- 
ferentiation, miR-124 is downregulated during 
cardiomyocyte differentiation of BM-derived MSCs via 
targeting STAT3, which is essential in the induction of 
transdifferentiation (44). On the other hand, miR-138 
(10) and miR-22 (45) are involved in MSC osteogenesis, 
and miR-30 (46) and miR-204 (1 1) in adipogenesis. In this 
study, the aforementioned miR-22 is upregulated in both 
osteogenic and adipogenic differentiated BM-MSCs. miR- 
22, as well as other miRNAs identified in Supplementary 
Figure S3, may be involved in common mesenchymal dif- 
ferentiation pathways, rather than being lineage-specific 
differentiation regulators. miRNAs associated with 
lineage-specific differentiation were pinpointed in the 
present study (Figure 3A), and their roles in MSC differ- 
entiation awaited characterization. Among these 
miRNAs, miR-19a-3p and miR-19b-3p, which are 
encoded by the miR-17-92 cluster, are increased in BM- 
MSC-derived adipocytes. Overexpression of miR-17-92 
could accelerate adipocyte differentiation in mouse 
preadipocyte 3T3L1 (47). Our results showed that miR- 
19a-3p and miR-19b-3p may be the critical miRNAs 
regulating MSC adipogenesis. 

CXCL12, also known as stromal cell-derived factor- 1, 
is an a-chemokine that is also a chemoattractant 
for other stem cells such as hematopoietic progenitor 
cells. CXCL12 stimulates migration by rearranging the 
actin cytoskeleton, increasing focal adhesion and matrix 
metalloproteinases production in MSCs (48,49). CXCL12 
is significantly upregulated in MSCs exposed to tumor 
cell-conditioned medium compared with cells treated 
with control medium, which suggests that CXCL12 sig- 
naling is important to bring about MSC migration toward 
the tumor microenvironment (48,49). SIKE1 (suppressor 
of IKKe) is known to interact with IKKs and TANK- 
binding kinase 1 to bring about inhibition of virus-trig- 
gered and toll-like receptor 3 (TLR3)-triggered activation 
(50). This does not occur via IkB-oc degradation, but 
rather IKKe and TANK-binding kinase 1 regulate NF- 
kB activity through a non-canonical pathway that directly 
targets the p65 subunit of the NF-kB complex via phos- 
phorylation of serine 536 (38). In cancer cells, IKKe is 
involved in regulating constitutive NF-kB activity, but 
not cytokine-induced NF-kB activity (51). In the context 
of WJ-MSCs, we proved that miR-146a-5p is a new 
CXCL12 and SIKE1 inhibitor. Regulating miR-146a-5p 
expression in MSC may benefit engraftment of trans- 
planted MSC that is homing onto injured tissues. 

Injured tissues secrete CXCL12 to recruit MSC for re- 
generation (52,53). The CXCL12/CXCR4 axis mediates 
MSC migration through activating PI3K/AKT (54), 
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MAPK/ERK (48) and Jak/Stat (48) signaling pathways. 
We found that miR-146a-5p is a downstream of CXCL12 
(Figure 6C). Inhibiting CXCL12 expression resulted in 
miR-146a-5p upregulation. How CXCL12/CXCR4 sig- 
naling pathways regulate the expression of miR-146a-5p 
is still unexplored. Promoter of miR-146a gene contains 
potential binding sites of STAT, c-ETS, PUT and NF-kB 
(55-57). Expression of miR-146a-5p may be regulated by 
these transcriptional factors in response to CXCL12 
stimulation. Transformed B lymphoma cells possess 
increased AKT phosphorylation and oncomiR miR-21 ex- 
pression, as well as downregulated miR-146a levels (58), 
indicating the PI3K/AKT pathway may involve in the ac- 
tivation of miR-146a-5p. Another candidate circuit 
involved in miR-146a-5p regulation is the NF-kB 
pathway: there are two NF-kB binding sites within the 
miR-146a-5p promoter region, and in vitro reporter 
assays showed that miR-146a-5p is induced via a NF- 
KB-dependent manner (30). We found that knockdown 
NFKB1 resulted in downregulated miR-146a-5p in 
WJ-MSCs (data not shown). 

The uncoupled MSC motility and growth in vivo may 
enable MSCs to migrate to regenerative sites via circula- 
tion at a single cell level while avoiding getting stuck in 
small blood vessels, hence benefitting MSC-mediated 
tissue repairing. Among MSCs, migrated cells have been 
shown to consist of more cells in G0/G1 phase and less 
cells in S and G2/M phase compared with non-migrated 
cells (59). Activation of NF-kB activity by TLR ligands is 
able to increase MSC proliferation but decreases MSC 
basal motility (36). Comparing of migrating and non- 
migrating fetal BM-MSCs showed that the nuclear recep- 
tors, Nur77 and Nurrl, can increase migration of MSC, 
while decreasing proliferation (60). The relationships 
between cell cycle and cell motility have been investigated 
in other stem cells or cancer cells. In HSCs, inhibition of 
the cell cycle may enhance HSC transplantation efficiency 
in vivo (61). Transforming growth factor p (TGF-P) and 
bone morphogenetic protein 4, on the other hand, mediate 
growth arrest and increase in cellular migration and 
invasion, which may eventually benefit distant tumor me- 
tastasis (62,63). Our results may also be applied to cancer 
research, as miR-146a-5p is a known migration inhibitor 
in various types of cancers (35,64). 

In summary, this study provides comprehensive 
miRNome profiles of two different types of MSCs, as 
well as cells differentiated from BM-MSCs. We envision 
that our genomic data will help to create a roadmap to a 
better understanding of MSC sternness and the regulation 
of MSC differentiation. Because both types of MSCs have 
potential as sources of stem cells for clinical transplant- 
ation, our findings have implications in terms of the 
clinical applications of MSCs. A manipulation of the 
level of miR-146a-5p in transplanted MSCs may benefit 
cell-based therapy in the future. 
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